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Abstract 

We analyze consequences of the approximation u\ — /x^ -C (a 'BPS' limit) for B and 
D mesons. It is shown that neglecting perturbative effects many power corrections would 
vanish to all orders in 1/tuq, in particular those violating heavy flavor symmetry. Among 
them are corrections to B — > D formfactors. A number of relations receive corrections 
only to the second order expanding around the BPS limit to any order in l/mg, including 
both / + and /_ at zero recoil. This allows an accurate evaluation of T + for B^D Iv. Its 
perturbative renormalization is computed analytically in the required Wilsonian scheme, 
yielding the dominant 3% enhancement. 
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The heavy quark expansion is a powerful tool for treating strong interactions relying 
on the first principles of QCD. Its indispensable part is a general Wilsonian separation 
of effects originating at sufficiently small and larger distances, a general concept of the 
operator product expansion (OPE) [1, 2]. The resulting expansion is most informative 
for inclusive decay probabilities where the OPE in local heavy quark operators emerges 
[3, 4]. Comprehensive application of all elements of the heavy quark theory often allows 
to put nontrivial constraints on the nonperturbative parameters in individual transition 
amplitudes based on information from inclusive decays. 

An important role here is played by the heavy quark sum rules, including recently 
suggested spin sum rules [5] for the Small Velocity (SV) regime [6]. For example, they 
put a rigorous bound g 2 > | on the slope of the IW function and lead to a number of 
other inequalities in the heavy quark limit. The recently proposed generalizations to 
higher orders in velocity [7] resulted in remarkable relations (D'Orsay sum rules) and 
nonperturbative bounds for its higher derivatives. The constraining power of the whole 
set of the SV sum rules depends strongly on the actual size of the B meson expectation 
value (/% of the leading local heavy quark nonperturbative kinetic operator Q{iD) 2 Q. 
Its theoretical expectations used to be a controversial subject for years. The inequality 
between the expectation values of properly defined kinetic and chromomagnetic operators 
puts a nontrivial lower bound on /i 2 . 

Although strong dynamics at small and large distances is governed by the same QCD 
equations of motion, physics originating from below and above a GeV scale is quite dif- 
ferent. It has been pointed out [8] that experiment may implicate an interesting pattern 
for the nonperturbative domain in B and D mesons, by favoring (see, e.g. Ref. [9]) fi 2 in 
the lower part of the allowed domain, 1 only little exceeding chromomagnetic expectation 
value //g(lGeV) = 0.35±;q2 GeV 2 . In this case it is advantageous to analyze nonpertur- 
bative dynamics combining the heavy quark expansion with expanding around the point 
where /i 2 = /Iq would hold. 

This is not just an arbitrary point of a continuum in the parameter space, but a 
quite special limit where the heavy flavor ground state has to satisfy functional relation 
an\B) = 0, with a and ir being heavy quark spin and momentum, respectively. From this 
perspective it is reminiscent to the 'BPS '-saturated state like the lowest Landau level for 
an electron in magnetic field. 2 It is not clear how deep the analogy with BPS symmetry 
goes, for the Pauli Hamiltonian describes the leading-order power corrections rather 
than the static heavy quark Hamiltonian itself 

Hoc = flight - / d 3 x Q + A Q (x) (1) 

shaping the ground state. The functional relation would not be generally respected by 
short- distance perturbative exchanges. It can only be viewed as an approximate property 

1 Certain theoretical loopholes in the framework applied in individual analyses have been identified 
[10, 11, 12] which are expected to be eliminated in the forthcoming experimental results. 

2 The analogy to this problem in quantum mechanics and possible relation to BPS symmetry were first 
noted by M. Voloshin (1999, private communication), whom the suggested name of the limit ascends to. 
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of nonperturbative dynamics at energies below or about 1 GeV scale, and applied to the 
ground state. 

A practical consequence of the proximity to the BPS regime is small room the heavy 
quark sum rules leave for possible values of the IW slope 

- £ g\l GeV) £ 0.95 at GeV) £ 0.45 GeV 2 . (2) 

A later UKQCD lattice evaluation [13] g 2 = 0-83-.\i-x>i fit well that prediction, although 
precision remains insufficient. This and some other consequences were noted in Ref. [8] (see 
also [10, 12]). In the present Letter we examine this remarkable limit in more detail. We 
will largely abstract from the perturbative corrections, in particular those renormalizing 
power corrections to the heavy quark Lagrangian. To elucidate underlying physics we 
often use quantum-mechanical language; correspondence with second-quantized notations 
in field theory is given, for instance in Ref. [14]. The necessary introduction into heavy 
quark expansion technique and used notations can be found there, as well as in reviews 
[15, 16]. One of the practical applications of the 'BPS expansion' is an accurate model- 
independent determination of \V c b\ from the B — > D Iv rate near zero recoil if enough 
statistics is available for this mode. Decays B^Dtu t sensitive to possible Higgs effects 
add motivation for a precision control. Power corrections to the decay amplitudes in this 
kinematics are obtained through order l/rng. 

1 'BPS' relations for B and D mesons 

The starting consequence of the BPS limit is equality of the spin-singlet and spin-nonsinglet 
expectation values appearing in the SV sum rules for asymptotically heavy quarks: 

Q 2 = \, A = 2S, pls = -pl- (3) 

Similar relations hold for nonlocal correlators of the 1/rriQ terms in the heavy quark 
Lagrangian: 

pIg = -2/4n p\ + pIg = ~(pL + pI), (4) 

a series extending to higher-order correlators. 

The BPS limit actually generalizes the heavy flavor symmetry to all orders in I/thq. 
First, all terms oc l/rng in the 1/mg expansion of the effective Hamiltonian annihilate 
the ground state. Then all corresponding T-products likewise vanish in it, and no power 
corrections to the heavy quark relation M P = mq + A appear, so that 

M B - M D = m b - m c (5) 

holds to all order in 1/mg. It is also important that the Foldy-Wouthuysen transfor- 
mation acts trivially (is unity) on the ground state, therefore the proper nonrelativistic 
wavefunction coincides with the upper component of the full Dirac bispinor in the meson 
restframe. 
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The above facts follow from the observation that for the 'BPS' wavefunction the lower 
component of the bispinor does not appear even when power correction are included. To 
prove this, we recall the full QCD equation of motion for the quark field in the nonrela- 
tivistic notations. If the static solution ip for mg — > oo, it (p = additionally satisfies 
the 'BPS' constraint (<77r)</? = 0, the bispinor 



Q(x,t) = e-^ t i (6) 




solves the Dirac equation 

( i + ™q1o)Q — m Q Q ? n^ = iD^- rriQV^ , (7) 
and the corresponding wavefunction 

si, (* ( T n _ ( K(xQ, {flight}) a = l,2\ , , 

is a formal eigenstate with E = uiq + A of the finite-mg Hamiltonian including light 
degrees of freedom - their wavefunction simply does not depend on ttiq. 

Absence of nontrivial Foldy-Wouthuysen corrections in the ground state is also trans- 
parent. Since time derivative in the QCD Lagrangian for fermions enters linearly with 
the gauge- field-independent coefficient, the only source of the transformation to nonrel- 
ativistic wavefunction is eliminating the lower component of the bispinor. If the latter 
does not appear, Eq. (6), no transformation is required. This corresponds to the fact that 
Foldy-Wouthuysen transformation Sfw (considered in the space of the upper components, 
in the restframe) can be viewed as 

S ™-{-^) " 1+ 8^+ 8^ +- ( 9 ) 

where the inverse square root of the projection operator is understood as the local OPE 

expansion in 1/mg of the forward matrix element Q 2 Q in terms of the upper 

components. As shown in Ref. [17], this expansion generates the whole nonrelativistic 
expansion for Dirac Hamiltonian, so triviality of Foldy-Wouthuysen transformation is 
intimately related to absence of corrections to the meson mass. 

Vanishing of corrections to heavy flavor symmetry in I/thq expansion applies also to 
the transition amplitudes between the ground-state mesons with different heavy quarks, 
like in B — > D decays. Below we take a closer look at the B — > D amplitude assuming 
both b and c are heavy enough to meaningfully apply the expansion. Since axial current 
does not contribute here, we focus on cry^fe-induced amplitudes; the results apply to other 
allowed Lorentz structures as well (say, scalar relevant for charge Higgs contributions). 

In general, the B^D amplitude is described by two vector formfactors 

(D(p 2 )\c^b\B( Pl )) = f + (f)( Pl+P2 ) u + f_(f)( Pl - P2 ) u . (10) 
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At zero recoil, q = a single amplitude, viz. u = remains: Jo = (M B + M B .)f + (0) + (M B — 
M D )f_(0). In the heavy quark limit one has 

M B + M D M B -M D 
f +<y0 ^ ~ 2VM b Md ~ ~ 2\J M B M D ' ^ 

It has been noted [10] that in the BPS limit all power corrections to J vanish. A stronger 
statement applies: the /+ formfactor determining all decays amplitudes with massless 
leptons, at zero recoil keeps its asymptotic heavy quark limit value 

„ . . M B + M D , . 

/+(0) = 27mm ■ (12) 

Moreover, at arbitrary momentum transfer the heavy quark relation between / + and /_ 
remains valid in higher orders in I/tuq: 

'-<«"> "^w^- < 13 > 

Furthermore, dynamic power corrections in B^D amplitude vanish: 

2 M B + M D / M%+M 2 n _ q 2 \ 
J+[q ' \ 2MbMd ) ' 

where C( v a v u) i s the normalized formfactor in the infinite mass limit (IW function). 
Therefore, the B^D differential decay rate would more or less directly measure the IW 
function. Below we show a way to derive these BPS consequences. 

To relate /+ and /_ in Eq. (13) we can consider two independent amplitudes, C70& and 
id^J^, applying the QCD equation of motion we have 

id^Jn = (m b -m c ) cb , (D\idMB) = (M 2 B -M 2 D )f + + g 2 /- • (15) 

Both C70& and cb currents do not mix upper and lower components for bispinors. There- 
fore, if the initial B meson is at rest and its lower component vanishes, the two currents 
coincide up to the factor m b — m c leading to 

(m b -m c ) i(M B + E D )f + + (M B -E D )f_\ = (M 2 B -M 2 D ) /++ (M 2 B +M 2 D -2M B E D ) /_ . (16) 

Replacing m b — m c by M B — M D we arrive at the stated relation. Hence, it is a 
consequence of vanishing lower components and of absent power corrections to the meson 
masses in the BPS limit. 

Non-renormalization of the (exact) zero-recoil amplitude can be readily seen, for in- 
stance, in the heavy quark sum rule for Jo extended to arbitrary order in I/itlq: the r.h.s. 
of the sum rule would not get corrections, 



+ VlFl 2 = l-^£p i_) 2 _ p 3 d+pIs (J- + J-) (l L] 2 _ (17) 

t-~t 1 11 4 \m c m b J 4 \m c m b J \m c m b J — ■> \ J 
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and inelastic transition amplitudes likewise vanish order by order in the BPS limit. The 
simplest way to understand this, however is using its quantum mechanical meaning re- 
vealed in Ref. [18]. 3 We have seen that no corrections to the quantum mechanical wave- 
function appear, lower components are absent and there is no explicit corrections in the 
C70& current - then no room for power corrections remains. 

Absence of corrections to the heavy quark limit relation (14) can be understood in the 
following way, looking once again at the timelike component of the current. Assuming 
B at rest and D moving with fixed velocity, we observe that the amplitude does not 
depend on m b - lower component of the b field is absent, while the upper one remains 
mass-independent. Therefore, it can only be a function of m c and can be computed at 
m^^oo. Changing the roles of b and c we would find it rather can only be a function of 
mj. Therefore it has to be free from any power corrections, at a given mesons velocity. 

Thus, the heavy flavor symmetry for the ground state would extend to any order in 
1/tuq expansion in the strict BPS limit. 

What if rn Q -> ? 

Absence of all power corrections to heavy quark relations even in a very special regime 
may sound paradoxical. We definitely know that when mq — > usual B and D mesons 
would rather evolve to become counterparts of the Goldstone particles like ir and K. The 
mass relation for them would be clearly different, for instance. 

There is no contradiction between the two regimes even if we abstract from the fact 
that the BPS limit cannot be exact and the corrections to it analyzed in the 1/mg 
expansion would explode when mq descends below a certain hadronic mass. The scale 
where the BPS-protected heavy quark relations get completely destroyed may even not 
decrease but remain stable when approaching the true BPS limit. For power expansion 
cannot have a finite radius of convergence being only asymptotic. Even if all power terms 
would vanish, there remain exponential terms scaling like 

e ^hadr (18) 

which modify the heavy quark relations. 

Presence of such effects is a general feature of expansions in field theory. In the 
heavy quark expansion they can explicitly appear due to unlimited spectral density of the 
operator 7To (heavy quark Hamiltonian), with the support stretching below — Ivtiq. Heavy 
quark expansion excludes extra heavy quark degrees of freedom expanding l/(2mQ + 7To) 
in l/2mQ. The expansion is convergent only if 1 7r 1 < 2rriQ in the operator sense, which 
seems impossible for a state well localized in space. Therefore, we should rather expect 
that violation of the BPS relations becomes of order unity below a certain hadronic scale 
regardless of proximity to the BPS regime for sufficiently heavy quarks. 

There are transparent mechanisms for such a change in the regime. An obvious candi- 
date is spontaneous breaking of the chiral symmetry through nontrivial condensate (QQ) 

3 Thc 1/jBq term for the sum rule for B^D* was derived in Ref. [15]. The simplest way, in fact uses 
this QM interpretation and explicit form of Foldy-Wouthuysen transformation, Eq. (9). 
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Figure 1: 'Level crossing' would completely change mQ-behavior below a certain hadronic mass 
even at arbitrary weak "BPS perturbation" 

of the QQ-p&irs in the physical vacuum when the quark becomes light. The basic as- 
sumption behind heavy quark symmetry, in contrast is separate conservation of numbers 
of Q and Q. Even if a formal solution of equations of motion exists, it may not be the 
correct ground state over the true vacuum. 

Another related mechanism is usual quantum mechanical level crossing. If analytically 
continued eigenvalues for two states happen to cross at some value of the heavy quark 
mass, the "repulsion" occurs, and behavior of the true ground state below and above the 
crossing mass becomes totally different, while the effect of any perturbation mixing the 
states becomes of order unity regardless of its formal strength, see Fig. 1. 



2 Expanding around the BPS regime 

The BPS limit cannot be exact in actual QCD, rather an approximate property of strong 
dynamics in the nonperturbative domain of momenta below a GeV scale. Therefore, like 
with the conventional I/tuq expansion it is important to determine the scale of violation 
of its particular predictions. If a concrete low-order in 1/mg correction is found in terms 
of the heavy quark operators, it is not difficult to see its BPS scaling. We, however, need 
a deeper classification which would hold to all orders in 1/rriQ. 

The practical expansion parameter to quantify deviations from the BPS limit, the 
norm of the state obtained by acting an on the ground state 



ll^is) II = V/4-/4, (is) 

has dimension of mass. A similar dimensionless parameter is 



P = ||7roV7?) 15)11 EE A /3 £? 2 - 



E 



|,W|2 
\ T l/2\ 



(20) 



(/T 1 in quantum mechanical notations is simply — — here). Numerically (3 is not too 
small, similar in size to generic l/m c expansion parameter in conventional l/m c series. 
Relations violated to order (3 may in practice be more of a qualitative nature, while 
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^2 ^ go can provide enough suppression. Moreover, we can count together powers of 
l/m c and powers of (3 to judge the real quality of a particular heavy quark relation. 

BPS relations (13), (14) for the decay amplitude at arbitrary recoil indeed receive 
corrections oc/3 1 , likewise equality of p\ G and — 2p^. 7r . Other relations between the heavy 
quark parameters mentioned in Sect. 1 are accurate up to terms f3 2 . In particular, 

• Mb — Mr, = mb — m c and Mr, = m c + A. 

• Zero recoil matrix element (D\cry b\B) is unity up to f3 2 . 

• Experimentally measured B — > D formfactor / + near zero recoil receives only second- 
order corrections in (3 to all orders in I/uiq: 

This is an analogue of the Ademollo-Gatto theorem for the BPS expansion. The similar 
statement clearly applies to /_ as well. 

Absence of /3 1 corrections to the meson mass is easily seen in usual perturbation theory: 

M B = (B\H tot \B) . (22) 

The terms linear in (3 would reside either in \B) or in 7i to t- The former vanishes since 
to the leading, (3° order \B) remains the same BPS state at finite m b . Absence of the 
linear in (3 terms from 7i tot is also understood - all local operators there include left-most 
an acting on b and right-most an acting on b since originate from excluding the lower 
bispinor components; this is transparent in the approach described in Ref. [17] . 4 

Absence of linear non-BPS effects in zero-recoil (D\cry b\B) is also transparent in BPS 
perturbation theory. Corrections oc/9 can originate either from charm or from beauty 
wavefunction, but not from both simultaneously, 

Spi(D\cy b\B) = (6(,i* D \*° B ) + (* D \6p* B ) = , (23) 

where we have used that the current acts on the unperturbed BPS state as a unit operator 
except changing flavor. Alternatively, this follows from the zero-recoil vector sum rule as 
described in the preceding section - power corrections to the sum rule appear only to 
order [3 2 , and excitation transition amplitudes appear to order (3. 

Vanishing of 0(/3) corrections in Eq. (21) for /+(0) is least obvious. Since this holds 
for the combination of /+ and /_ describing J , it suffices to show this for an alternative 
combination, say taking the spacelike current cryb in the case where one of the mesons is 
at rest and retaining only linear in velocity terms. Once again we make use of the fact 
that order-/? 1 effects could originate from perturbation in either B or D wavef unctions, 
but not both simultaneously. To simplify algebra we can assume that B meson is at rest 
if BPS corrections appear in B, and go to the D restframe if a deviation from the BPS 
state occurs in charm sector. 



4 In particular, this follows from its Eq. (10) once the scalar expectation value (B\bb\B) is 1 + 0{f} 2 ). 
The energy-momentum tensor V does not depend explicitly on Wf, and plays the role of the rest-frame 
Hamiltonian of light degrees of freedom. 
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The next helpful observation is that if, say the beauty sector wavefunction enters to 
the leading order in f3, it can be assumed to be in the more familiar mj, — > oo limit. 
Projecting the matrix element on v we then get the spinor structure 

"-^(sy^ ?)($?). <*> 

where ip^ and refer to -D at rest and include all mass corrections. 

The term (<// c )|(<7n)^|</?Q^) yields just the heavy quark limit result, Eq. (11) if ip^ 

(c) 

coincided with the asymptotic wavefunction (p . It generally does not, but the first-order 
corrections to the overlaps like {^Po^l^Po^} always vanish. This applies to BPS perturbations 
since (p^ remains normalized to unity up to terms oc/3 2 to any order in l/m c . 

Another source of power corrections is the product (x^\(^n)\ l Po'\v)). Now we use 
that (fo\v) is the asymptotic heavy quark wavefunction which velocity dependence to 
the first order is given by 

\<P { o\v)) = \<p$\o)) - -1^= nv\^\0)) = \<p$\o)) + n l nv\^\0)) , (25) 

/Too — A 

where Tioo is the static heavy quark limit Hamiltonian, and the second form uses the 
analogue of the second-quantized notations. Since heavy quark spin a commutes with 
Hoo this piece reads 

f(D\c^\an)b\B (v)) (26) 

when calculated for spinless meson states. Therefore it vanishes for the BPS beauty state 
regardless of the charm wavefunction. 

The above reasoning essentially relies on linear in v approximation. To order v 2 , for 
instance, the operator n 2 appears along with (<?7r) 2 , which does not vanish when acts 
on the BPS state. Therefore, while having managed to leave the point of zero recoil, we 
cannot extend the Ademollo-Gatto theorem to, say the slope of the formfactor. 



3 Application to B — ► D lis near zero recoil 

Since /_ formfactor does not contribute to any decay amplitude with massless leptons, 
the amplitude even near zero recoil depends on the space-like current which suffers from 
linear l/m c effects in the 1/mg expansion; conventionally this is viewed as the serious 
theoretical drawback of such a decay mode. The BPS expansion turns out more robust 
in this respect protecting the heavy quark relation up to the second order. This allows 
an accurate estimate of the B-^D rate near zero recoil in terms of | V^f, | 2 . 

Similar to heavy quark symmetry itself, the BPS limit is affected by short-distance 
perturbative effects. Accounting for the latter should therefore be done accordingly. The 
Wilsonian procedure with the explicit 'hard' cutoff \i is most suitable here simply elimi- 
nating the low-momentum domain. The principal perturbative corrections are just short- 
distance renormalization of the bare cry u b current itself. The technique for such Wilso- 
nian calculations has been elaborated and applied to zero-recoil B — > D* amplitude, see 
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Ref. [16]. The method specifically corresponds to the adopted renormalization scheme 
for nonperturbative operators with the upper cutoff \x on energy [8]. The case of B^D 
assuming non- vanishing velocity requires only technical modifications. 

The one-loop result with arbitrary cutoff scale /x is obtained in an analytic form: 



2a s 



3vr 



2{mi —m 



J 4 m c u b -m h Lo c _ _. h 



3ii 2 mfe— m c 



-| J£L_J±k. o™ ^ J L 1 1 



3 m b — m c m c + nib mb~m c 6 m c +mb \m c 



3 mb — m c 



+ — 13 



m. 



— '-r+-f +? / •> 2 > 



3m c + mb \uj^ u)% J 6 \m c rribJ 3m c mbmf ) —m 2 



where 



\im 2 c + /r. u b = 

Fig. 2 shows it numerically assuming m c =1.2GeV, mf, = 4.6GeV and a s = 0.3. 



(27) 
(28) 




Figure 2: Short-distance renormalization for B —> D amplitude at m c = 1.2GeV, mf, = 4.6GeV 
and a s = 0.3 as a function of the separation scale \i. 

The first-order 1/mq correction is readily read off the current in Eq. (24), cf. Eq. (26) : 5 

M B -M D 



* = T^HT >♦«» = + H (-"-J WTW ~ ° , 

M B + M D \ 2 J \m c m b J M B + M D \mf 



(29) 



it is positive, explicitly suppressed by two powers of (5 in accord with the theorem in the 
previous section, and is small numerically. One finds that linear p/rriQ dependence of 
one-loop £y(/i) cancels against the one-loop /i-dependence of f — S, as expected. 

Although 4 — S normalized at 0.8 to 1 GeV has not been directly measured, a good 
estimate for it would be the actual upper bound [16] 

-(B\l(3H)(-^ l )(39)b\B) < 



2 3 2M B 



2 _3\ /4 

4 / Q 



(30) 



3 A connection between the 1 /tuq correction and inclusive transitions was first noted in Ref. [19] . 
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The IW slope q 2 has not been well determined experimentally yet. An alternative estimate 
would use the scale of an average excitation energy e of the P-wave ^-states saturating 
the sum rules: _ 

— — E « (31) 
2 3e K ' 

with e 500 to 700 MeV. Depending on the precise value of /z 2 (1 GeV) the 1/m correction 

to JF + emerges at 1% level. 

Since the expansion in l/m c is involved, higher-order corrections a priori could be 
significant. This is illustrated by the perturbative contribution where linear in \x approx- 
imation works well only for \x ^ 400 MeV, while already at /i ~ 800 MeV second- and 
third-order terms are of the same size. The BPS expansion ensures, however that the 
overall suppression is carried on to all higher-order nonperturbative effects as well. 

The second-order power correction to JF + has three pieces, each manifestly of the 
second order in BPS. It also follows directly from Eq. (25), now explicitly expanding in 
l/m c : 

5 l/m ,F + =(—-—) {V x + V 2 + V 3 } (32) 

The first one, V\ is the same as in zero- recoil (D\cr/o\B) and is well understood [14]. It is 
negative consisting of local correction and nonlocal "overlap deficit" : 

V, = ■ (l + xZi) , (ASo) xZi = ^(B\b(<H0 2 tf(**) 2 b\B) 

= [i\x \d*x^(B\iTWair) 2 b(x),b(air) 2 b(0)}\By>0. (33) 
J AM. n 



IB 

The second piece is local and positive, V 2 = ^ ^ G . The last piece can be viewed as the 

1/mq correction to the value of ^ — £ in a finite- mass meson and is described by the 
expectation value of the nonlocal correlator 

~ 6Vs = ^( B \HW 2 *o\tt)*o\**)b\By . (34) 

It is plausible that the correlator is likewise positive decreasing the effective value of f— X, 
although strictly speaking the sign is not fixed. Its scale can be roughly estimated as that 
for the correlator Q(an)nQ 1 (a7v) 2 7TQ 1 (o : 7v)Q and taking for the latter &f {\ — ^) as an 
educated dimensional guess. Since local pieces nearly cancel, S 1 / m 2j 7 + is dominated by 
two nonlocal correlators. 

Collecting all terms through the second order in I/tuq we arrive at the estimate 

5 P owcr^+ £ 0.01 at /i 2 £ 0.43 GeV 2 , (35) 

where following Refs. [14, 20] we assume xZ-l~^-^ ^ 0-5- The literal prediction depends 
moderately on the actual kinetic expectation value, 

f£(l GeV) -0.43 GeV 2 , , 

F + ~ 1.04 + 0.13 ^ ; (36) 
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the uncertainty, however would soon go out of control at this level of precision for 
/4(lGeV) exceeding 0.45 GeV 2 . 

The above estimates suggest that the nonperturbative corrections to T+ are really 
tiny, and it can be accurately evaluated. It should be appreciated, however that practical 
implementations of the heavy quark expansion leave out possible exponential terms like 

<5ex P F+ oc ( e "hadr - e \ (37) 

which are routinely ignored, but may be essential at a percent level when relying on heavy 
charm expansion [21]. At the moment we cannot say much about them; yet they probably 
put the actual limit on the accuracy of the theoretical predictions for JF + if /^(i GeV) is 
confirmed to be below 0.45 GeV 2 . 



4 Conclusions 

We have analyzed the structure of the pure nonperturbative corrections to B and D 
mesons near the 'BPS' regime which would apply if nonperturbative physics is largely 
confined below 1 GeV scale and yield /z 2 close to /Iq. A number of relations are shown 
to hold extending the heavy flavor (but not spin!) symmetry to all orders in l/mg in 
the BPS limit, valid however only for the ground state pseudoscalar heavy flavor mesons. 
Some of the important BPS relations get corrections only to the second order in the 
BPS expansion regardless of the order in 1/wiq, among them are two B^D zero recoil 
formfactors. This is the analogue of the Ademollo-Gatto theorem for the BPS expansion. 

The practical utility of such an expansion strongly depends on the actual size of 
/x 2 (lGeV). If - as has been reported by most experimental analyses - it centers at or 
below 0.4, or even up to 0.45 GeV 2 , it is rather powerful. At larger values of /x 2 the 
predictability weakens; in this case, however any application relying on the heavy quark 
expansion for charm for non-BPS-protected relations would become questionable. 

Dynamic heavy quark expansion has enlightened us that the scale of nonperturbative 
effects for heavy quarks in actual QCD is quite significant, fi np > — 700 MeV at 
least. This contrasted early ideas that this scale is like a 'constituent' light quark mass 
^const ~ 250 MeV inherited from naive nonrelativistic quark models. Yet it turns out that 
the latter small scale can sometimes re-emerge in theory - where the 'BPS '-protected 
properties are considered, with — [x 2 G <^ 300 MeV. A better theoretical understanding 
of the dynamic origin of such a hierarchy is highly desirable. 

An interesting place to test the BPS predictions is the ratio /_(0)// + (0) fixed in the 
heavy quark limit, 

/-(0) = M B -M D 

/+(0) M B + M D ' 1 1 

Since both /+(0) and /_(0) undergo power corrections only to the second order in the BPS 
expansion, they must be small and accurately evaluable, Sect. 3, thereby offering a probe 
for possible exponential effects. The suppressed formfactor /_ can be measured in the 
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B — > D tv t decays in future high-statistics experiments. This mode also probes possible 
effect of charge Higgs exchanges, which requires a precision understanding of the SM 
amplitude. The BPS expansion provides support for such a treatment; the electroweak 
corrections are to be properly incorporated at this level. 

The heavy quark expansion together with heavy quark sum rules allows an accurate 
prediction for the B — > D amplitudes near zero recoil in terms of measured observables. 
Our estimate (the electroweak effects, in particular the factor of 1.007 from the univer- 
sal short-distance renormalization are not included here) depends on /^.(lGeV), but for 
moderate values the power corrections appear at a percent level: 

T + = 1.04 ± 0.01 powcr ± 0.01 pcrt + 5 exp , (39) 

with perturbative corrections from momenta above 1 GeV contributing the dominant piece 
of 3%. It can be further refined. The corrections are significantly smaller and more definite 
compared to the 'gold plated' B^D* decay mode. 

This assessment differs from the existing estimates [22] although, in principle is com- 
patible with the predictions within their respective large error bars. The rationale is 
readily seen equating and counting together orders in conventional l/m c and in the 'BPS' 
expansion: Eq. (29) gives power corrections through the third order, with estimates in- 
cluded for the fourth-order effects. In contrast, no BPS-backup exists for B^D*, the 
corrections are significant in the BPS regime [8], and they are uncertain already to the 
leading order l/m^. 
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